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Abstract 



The evolution of non-potential perturbations to a current-free magnetic X- 
point configuration is studied, taking into account electron inertial effects as 
well as resistivity. Electron inertia is shown to have a negligible effect on the 
evolution of the system whenever the collisionless skin depth is less than the 
resistive scale length. Non-potential magnetic field energy in this resistive MHD 
limit initially reaches equipartition with flow energy, in accordance with ideal 
MHD, and is then dissipated extremely rapidly, on an Alfvenic timescale that 
is essentially independent of Lundquist number. In agreement with resistive 
MHD results obtained by previous authors, the magnetic field energy and 
kinetic energy arc then observed to decay on a longer timescale and exhibit 
oscillatory behavior, reflecting the existence of discrete normal modes with finite 
real frequency. When the collisionless skin depth exceeds the resistive scale 
length, the system again evolves initially according to ideal MHD. At the end 
of this ideal phase, the field energy decays typically on an Alfvenic timescale, 
while the kinetic energy (which is equally partitioned between ions and electrons 
in this case) is dissipated on the electron collision timescale. The oscillatory 
decay in the energy observed in the resistive case is absent, but short wavelength 
structures appear in the field and velocity profiles, suggesting the possibility of 
particle acceleration in oppositely-directed current channels. The model provides 
a possible framework for interpreting observations of energy release and particle 
acceleration on timescales down to less than a second in the impulsive phase of 
solar flares. 



1 Introduction 



Renewed interest in the long-standing problem of energy release in solar flares 
has been kindled by the launch of the Ramaty High Energy Solar Spectroscopic 
Imager (RHESSI) spacecraft (Lin et al. 2002) and, very recently, by observations 
of exceptionally large flares (Bowler 2003). Resistive magnetohydrodynamics (MHD) 
provides the framework for most theoretical studies of this problem (for a recent 
review, see Priest & Forbes 2000). Craig & McClymont (1991, 1993) and Craig & 
Watson (1992), for example, have invoked the resistive MHD relaxation of a two- 
dimensional magnetic X-point as a paradigm for energy release in flares. The magnetic 
equilibrium configuration in this analysis was potential (and thus stable), and the 
system was assumed to have a circular boundary centered on the X-point. Craig & 
McClymont (1991, 1993) investigated the stable mode spectrum of such a configuration. 
A countably infinite set of discrete normal modes was identified, with azimuthally 
symmetric modes corresponding to topological reconnection. Both the real frequency 
and the damping rate were found to increase with the number of radial nodes. It 
was inferred that the spectrum after a sufficiently long time would be dominated by 
the mode with the longest radial wavelength, i.e. the lowest damping rate, regardless 
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of the initial configuration. Tliis damping rate was sliown analytically to be given 
approximately by 

^"''^ - 2Ro{]nSy' 

where 5" is the Lundquist number at the radial boundary of the system R — Rq and 
Cao is the Alfven speed at that boundary, computed using the equilibrium magnetic 
field. The logarithmic dependence of the dissipation rate on S is reminiscent of that 
predicted by the Petschek "fast" rcconncction model (Petschek 1964), although it 
should be noted that Petschek used a steady state approach while the treatment of 
Craig & McClymont was time-dependent. The latter authors noted that, for typical 
solar parameters, equation (1) yields an energy release timescale of the order of several 
minutes to an hour, This is sufficiently rapid to explain thermal energy release in the 
gradual phase of a flare, but too slow to account for the impulsive phase: as noted 
by Tandberg-Hanssen and Emslie (1989), the thick target interpretation of flare hard 
X-ray bursts requires a field of the order of 0.02T to be completely annihilated in a 
volume of around lO^^m^ every second. 

Recognizing that resistive MHD may not provide a full description, several authors 
have investigated the possible role of non-MHD effects in fiare energy release. Craig & 
Watson (2003), for example, have identified exact reconnection solutions that include 
Hall current and electron inertial effects. Such effects have been studied extensively 
by researchers in the magnetospheric community: Birn et al. (2001) have summarised 
the results of a recent project to compare computational models of reconnection in 
Harris-type geometry ranging from resistive MHD to particle-in-cell. In a recent 
paper McClements & Thyagaraja (2004) generalized the spectral analysis of Craig 
& McClymont (1991, 1993) to include electron inertia. A key motivation for doing so 
in the laboratory context is that reconnection events in magnetic fusion experiments, 
such as those responsible for "sawtooth" oscillations in tokamaks, cannot be explained 
using MHD alone. As noted by Wesson (1991), the observed relaxation time in a 
sawtooth crash can be much shorter than the timescale implied by purely resistive 
models. McClements & Thyagaraja used the following form of Ohm's law: 

E + vxB = r/j + ^^, (2) 
ne^ at 

where rj is resistivity, n is particle density, and me, e denote the electron mass and 
charge. McClements & Thyagaraja demonstrated that the magnetic X-point eigenmode 
spectrum in this case has a continuous component in addition to the discrete spectrum 
studied by Craig & McClymont. This continuum replaces the ideal MHD Alfven 
continuum which is present when both resistivity and electron inertia are excluded. 
All of the finite frequency continuum modes have the same intrinsic damping rate, 
namely 

_ Cap 
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where 5e = c/ (UpeRo), c being the speed of Ught and the electron plasma frequency, 
is the coUisionless skin depth normalized to the system size. Because 6^ is typically 
very small, the continuum mode damping rate given by equation (3) can be high in 
absolute terms. In fact, if the resistivity is determined by electron-ion collisions, 7cont 
is essentially equal to the electron collision frequency 



where rii^ is particle density in units of lO^^m""^ and Tq is electron temperature 
in units of lO^K. Thus, for typical solar coronal parameters and Spitzer resistivity, 
the energy dissipation timescale implied by equation (3) is a fraction of a second. 
Statistically significant fluctuations have been observed in flare hard X-ray emission 
on such timescales (e.g. Kiphnger et al. 1983): this suggests that non-MHD effects 
may be significant in the context of flare energy release. However, in order to compute 
dissipation timescales reliably, it is appropriate to solve an initial value problem rather 
than using an eigenvalue approach. Craig & Watson (1992) performed such an analysis 
using a resistive MHD model: in this paper we adopt a similar approach, but with 
the electron inertial term in Ohm's law taken into account. Speciflcally, we solve a 
set of linearized equations describing a perturbed magnetic X-point using parameters 
that range from the resistive MHD limit considered by Craig & Watson (1992) to 
the coUisionless limit. Our approach differs from the two-fluid analyses carried out 
for example by Biskamp, Schwarz, & Drake (1997) and Ramos, Porcelli & Verastegui 
(2002) in that we focus on a single, dominant non-MHD effect, namely electron inertia. 
Compared to resistive MHD, the problem is then characterised by only one additional 
parameter: the coUisionless skin depth. Using this model, we address the key issue 
of how rapidly the energy in a non-potential magnetic fleld perturbation is dissipated 
or converted to kinetic energy, making an exact comparison with the resistive MHD 
scenario investigated by Craig & Watson (1992). 

This paper is structured as follows. After formulating the linearized initial value 
problem in §2, we proceed to solve it analytically and numerically for various parameter 
regimes in §3. The possible relevance of these solutions to energy release in solar flares 
is discussed in §4. 

2 Formulation of Initial Value Problem 

We consider inviscid incompressible perturbations of a two-dimensional current-free 
magnetic X-point in the low plasma beta limit. With Ohm's law given by equation 
(2), the induction and momentum equations can be written in the form (McClements 
& Thyagaraja 2004) 



i/e ^ 30 
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^ + {^r■V)v^- — {V'^^;)VA (6) 
at jiop 

where is a magnetic flux function whose curl yields the magnetic fleld, p and v denote 
fluid density and velocity, and //q is vacuum permeability. Both the magnetic fleld and 
V are assumed to lie in the (x, y) plane of the X-point. We linearize these equations 
by neglecting terms of second order in v and putting il) — ipE + ip where 

■0 is assumed to be azimuthally symmetric, and |VV'| is assumed to be much smaller 
than IVt/'eI- In equation (7) Bq is the unperturbed magnetic fleld in the {x,y) plane 
at radius R = Rq and 9 denotes azimuthal angle. When linearized, equations (5) and 
(6) become 

|(^_|v^^)+(v.V)^, = ^V=^. (8) 
W = 

where po denotes the unperturbed density (assumed to be uniform). These equations 
can be combined to give a single equation for ip in which 6 does not appear except in 
the Laplacian operator (Craig & McClymont 1991; McClements & Thyagaraja 2004). 
In the linear approximation, the variation of ip with 6 can thus be separated from the 
variation with R and t. In this paper we concentrate on the azimuthally symmetric 
case in which ip is independent of 9. 

It is clear from equations (7) and (9) that the components of v are ^-dependent even 
when ip is azimuthally symmetric. Formally integrating equation (9) with respect to 
time and writing the result explicitly in terms of {R, 9, z) components, we obtain 

BR f ~ 

v = — ^ (cos 2^, -sin 2^,0) V'^i^dt = v{cos29, - sm29,0). (10) 

The function v{R,t), deflned in this way, satisfles the scalar momentum equation 

I = ^V^^, (11) 
dt PqPqRq 

and the radial and azimuthal velocity components are given by 

vr = vcos29, (12) 

V0 — —vsm29. (13) 

It should be noted that v in these expressions can be either positive or negative. The 
induction equation then becomes 

. !^ = (14) 
/to po 
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Following Craig & McClymont (1991), we make the system dimensionless by 
normalising R to i?o, v to caq = Bq/ {hqPq)^/'^ , t to Rq/caq and ip to BqRq. Dropping 
the tilde on we thus obtain 




) 



vr + —— r— — 



Sr dr \ dr 



1 d f dip 



) 



(15) 



dv d f dip 
dt dr \ dr 



(16) 



where S = iiqRoCao/v- We proceed to solve these equations subject to initial conditions 
on dijj/dr and v, and boundary conditions on dip /dr. Regularity requires that the 
perturbation to the azimuthal magnetic field —dip/dr vanish at r = 0; this is consistent 
with the linear approximation since the equilibrium field in the (x, y) plane has a null 
at the X-point. We also set dtp / dr = Q r = 1, for reasons that will be explained 
in the next section. Following Craig & Watson (1992), we consider both spatially 
extended initial field perturbations and localized ones. A simple obvious choice of 
initial condition for adopted throughout this paper, is = 0. 

3 Solution of Initial Value Problem 
3.1 Energy Evolution 

Our principal objective is to compute the evolution of field and kinetic energy associated 
with a perturbed X-point. In the general case, with finite collisionless skin depth 5e, the 
total energy 8 has three components: non-potential field energy £f, ion kinetic energy 
£]a and electron kinetic energy £ke. In terms of the dimensionless variables introduced 
in the previous section, these energy components can be represented by the integrals 



The identification of the integral in equation (17) as the non-potential field energy 
follows from the fact that the magnetic field perturbation is —dip/dr. Locally, there 
is a contribution to the field energy density proportional to {dips/ dr) {dip / dr)., but 
this gives a zero contribution to the total energy when integrated over 6 since dips/dr 
is proportional to cos2^ whereas dip/dr is azimuthally symmetric. Ions make the 
dominant contribution to the bulk fluid velocity v, and therefore it is appropriate to 




(17) 




(18) 




(19) 



6 



regard the integral in equation (18) as a measure of ion kinetic energy. On the other 
hand, most of the current is carried by electrons, and the current density (which is 
oriented in the z-direction) is proportional to —'V^ip: it follows from this that the 
integral in equation (19) represents electron kinetic energy. 

In the limit S ^ oo the total energy 

1 /•! 

2 Jo 

is conserved. To estabhsh this result, we evaluate the time derivative of £, using 
equations (15) and (16) to eliminate dt/j/dt and dv/dt, obtaining 

f^-i/'(vV)Vd. (21) 

where we have invoked the boundary conditions dijj/dr = at r = and r = 1. 
The right hand side of equation (21) is simply the rate of Ohmic heating. By setting 
dt/j/dr — Bit r — 1, we ensure that the Poynting flux is locally zero at all points 
on the boundary. We are not setting ^; = at r = 1, and therefore we are allowing 
the possibility of a local mass flow through the boundary. However, since the radial 
component of the velocity vector, like the equilibrium field, has a cos2^^ dependence 
(eq. [12]), the integrated mass flux through the boundary is always zero. Thus, after 
t — there is no net flow of energy into or out of the system, and the energy can only 
decrease, through Ohmic dissipation, with the rate of decrease determined purely by 
the internal dynamics of the system. 

It is clear from equation (21) that the system is conservative in the limit 5" — > oo, even 
when 5e is flnite: electron inertial effects are essentially reactive and thus do not give 
rise to dissipation. However, we will demonstrate that the field energy Sf can change 
rapidly when 6^ is finite and S ^ oo. The conservation of S in this limit provides a 
critical test of the numerical schemes used for solving the full equations. In Appendix 
A we obtain a relationship between field energy and electron kinetic energy that will 
prove to be useful for the interpretation of finite 5e solutions presented in §§3.4. 

3.2 Ideal MHD Solutions 

As first noted by Bulanov & Syrovatskii (1981), a simple analytical solution of 
the linearized equations exists in the ideal MHD limit {S — > oo, 5e — 0). This 
solution provides an additional test of the numerical schemes, and moreover is key 
to understanding certain effects observed in the non-ideal case. 

Following Craig & Watson (1992), we evaluate the ideal MHD solution explicitly for 
two alternative initial conditions. When Sg — 1/S — equations (15) and (16) can be 



dr 



+ v^ + 6l{\/^-tPf 



rdr, 



(20) 
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combined to give 



With the substitution u = Inr, equation (22) reduces to the one dimensional wave 
equation 



with general solution 



(23) 



^ = f(^u + t)+ g{u -t) = /(Inr + t) + c/(lnr - t), (24) 

the functions / and g being arbitrary. Employing the boundary conditions dijj/dr — 
at r = and r = 1 and initial conditions 

— = sm(7rr), (25) 

it is straightforward to show that / and g are such that the full solution is given by 

■0 = — 77- |cos(7rre*) + cos(7rre~*)| , t < — Inr (26a) 

■0 = jcos (^-e *^ + cos(7rre~*) I , i > - Inr (26b) 
Putting 5e = 1/ S = Q in equation (15), we obtain 

■u = ^ |e*sin(7rre*) - e~*sin(7rre"*)| , i < - Inr (27a) 

-Inr (27b) 




The solutions for f > — In r represent smooth continuations of those for t < — In r: there 
is no discontinuity at this point. Using these expressions, we can compute the time 
evolution of the ion kinetic and field energy. As shown in the left plot of Figure 1, the 
two components of the energy undergo a single oscillation before reaching equipartition 
in about three Alfven times. The field and velocity profiles continue to evolve after 
this time, however. As shown in the right hand frame of Figure 1, the field energy 
and kinetic energy become increasingly concentrated in the vicinity of the X-point. 
Eventually, the field and velocity gradients become sufficiently large that the neglect of 
resistive and electron inertial terms in equations (15) and (16) is no longer justified and 
the ideal MHD solution becomes invalid. The right hand frame of Figure 1 exemplifies 
the well-known tendency of magnetic X-points to focus and accrete electromagnetic 
and kinetic energy (see e.g. Craig & Watson 1992). 
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The trajectory of the inward-propagating wave described by equation (24) is given by 
dr/dt = — r: this arises from the Alfvenic character of the wave and the fact that 
the magnitude of the equihbrium magnetic field is proportional to r (cf. eq. [7]). 
McLaughlin & Hood (2004) have recently demonstrated an important consequence of 
this, namely that a disturbance initially consisting of a plane wave is refracted as it 
approaches the X-point in such a way that it becomes more azimuthally symmetric: the 
region of the wave front closest to the null propagates more slowly than neighboring 
regions. For this reason, it is particularly appropriate to consider the evolution of 
azimuthally symmetric perturbations. 

It is also instructive to compute the ideal MHD solution for the case of a more localized 
initial flux perturbation of the form 



exp 




(28) 



where Ait is a constant characterizing the width of the disturbance in In r space. Craig 
& Watson (1992) considered initial perturbations similar to this, but with a factor of 
Inr outside the exponential (so that ijj rather than dip/dr was set equal to zero at r = 1) 
and Am = 1. Profiles of this type are only strongly localized if Am ^ 1. Determining 
ip{r,t) for this case from equation (24), and using equation (15) with Se = 1/S = to 
compute V, we obtain 



dr 



2r 



(In r + t) exp 



Inr + t 
Au . 



+ (In r — t) exp 



Inr — t 
Au 



, (29a) 



V = — <^ (kir + t) exp 



Inr + 1 
Au 



(In r — t) exp 



Inr — t 
Am 



(29b) 



If Am ^ 1, dip/dr and v are generally negligible except for Inr ~ ±t. Since we are 
only considering t > and r < 1, i.e. Inr < 0, the dominant terms in equation (29) are 
those multiplied by Inr +t, representing an inward propagating wave. It is immediately 
apparent in this case that v ~ dip/dr, i.e. energy is equally partitioned between the 
non-potential component of the magnetic field and the bulk flow, as in the case of the 
global disturbance represented by equation (25). 

The energy components corresponding to this solution are plotted versus time for 
Am = 0.1 in the left hand frame of Figure 2. It is apparent that equipartition occurs 
much more rapidly than in the case of Figure 1. The right hand frame of Figure 2, 
showing the evolution of dip/dr, indicates as before rapid focusing of the field energy. 
Steepening of the profiles is noticeably more rapid than in the case of Fig. 1: we would 
accordingly expect non-ideal effects to occur at an earlier stage of the system evolution. 
The two solutions considered above illustrate the fact that energy equipartition is a 
natural asymptotic state of the ideal MHD system (cf. Craig & Watson 1992). 
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3.3 Resistive MHD Solutions 



In this subsection we present results obtained from numerical solutions of equations 
(15) and (16) for the resistive MHD case, Sg = 0. This scenario, with somewhat different 
initial conditions, was investigated extensively by Craig & Watson (1992): we consider 
it here in order to provide an exact benchmark for the finite 6e cases considered in §§3.4, 
and also to highlight a phase in the late evolution of S not discussed by those authors. 
Independently-written codes were used to solve equations (15) and (16) for 5e — and 
for finite Sg (§§3.4). The numerical methods used in these codes are described briefly 
in Appendix B. 

Our initial conditions are those defined by equations (25) and (28), together with 
V = 0, as before. In the case of the more localized field perturbation (eq. [28]), we 
take Au = 0.1. For the two types of initial field perturbation, we have determined 
the evolution of the system for a wide range of Lundquist numbers S, from 10 to 10^; 
Craig & Watson (1992) concentrated their investigation on the case of 5" = 10^. Our 
results, shown in Figures 3 and 4, concur at high S with those obtained by Craig 
& Watson. For S > 10^ there is an early phase in which the system evolution 
is well-described by ideal MHD. In Figure 3(h), for example {S = 10^), the time 
evolution of the energy components up to i ~ 7 is identical to that computed using the 
corresponding ideal solution (Figure 1). The ideal phase becomes progressively longer 
as S is increased: for very small S, on the other hand, there is no ideal phase and the 
energy decays immediately. Although the onset time of the drop in energy increases 
with S, the timcscalc of the decay itself is essentially independent of this parameter, 
and is comparable to or less than one Alfven time. The factor by which the energy falls 
at this point appears to depend critically on the initial field perturbation: this was also 
found by Craig & Watson. In the case of the more localized perturbation (Fig. 4), the 
energy is reduced in this phase by a factor of about ~ 2 x 10^. 

The rapid drop in energy is followed by a period of slower decay in which the field energy 
and kinetic energy exhibit oscillatory behavior. As noted by Craig & Watson (1992), 
these oscillations arise from the discrete spectrum of weakly damped eigenmodes 
studied by Craig & McClymont (1991). After a sufficiently long period, the oscillations 
disappear. This can be seen in Figures 3 and 4 for 5" < 10^, and in Figure 5 for 5" = 10^, 
10^. The oscillations persist for a time that increases as (InS*)^: this is consistent with 
the S dependence of the least damped discrete mode (eq. [1]). 

After the oscillatory phase, the energy continues to decay, but on a longer timescale, 
and the energy moreover is almost entirely kinetic {S]^i ^ S{). This late phase of non- 
oscillatory decay, which is apparent in Figures 3(a), 3(b), 4(a), 4(b) and 5, can be 
accounted for by the fact that the discrete modes identified by Craig & McClymont 
(1991) do not constitute a complete set: as demonstrated by McClements & Thyagaraja 
(2004), there is also a continuum of damped modes with zero real frequency. When Sg 
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and S are finite, equation (22) is replaced with the fourth order equation 





(30) 

ui' / o'l' ur \ ur i 
For %Ij = e~*'^*/(r), equation (30) reduces to 

.V/+(,.^-.^«-.|)A(4).o. (31) 

The coefficient of the highest order derivative in this equation vanishes for complex 
mode frequencies cu satisfying 

a;\^ + i| = r^ (32) 

McClements & Thyagaraja showed that this relation gives rise to two distinct continua: 
one with finite real frequency, the other with zero real frequency. The former, which 
we will discuss in the next section, only exists when Se 7^ 0. The zero real frequency 
continuum, on the other hand, exists whether 6e is finite or not. This can be seen by 
putting Lu = — i7, where 7 is taken to be real. Equation (32) then yields 
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1 ± (1 - Ar^S^S, 



(33) 



2S51 

In the limit S'^6l ^ 1, the negative root in this expression yields a damping rate 

-f = Sr\ (34) 

Thus, there is a purely damped continuum mode for each point in the solution domain. 

Since this extends to r = 0, the damping can be arbitrarily weak. In the resistive MHD 
case, the Fourier spectra of ip and v will include contributions from this continuum as 
well as the discrete modes identified by Craig & McClymont. Equation (34) indicates 
that non-potential perturbations to the X-point field will in general take an infinite time 
to be completely dissipated. The presence of the continuum only becomes apparent 
when the discrete modes have disappeared: since the discrete mode damping rate 
scales with 1/(^5*)^ (eq. [1]), the presence of the continuum becomes manifest at 
progressively later times as S is increased. 

As noted above, another feature of the late non-oscillatory decay phase is that the 
remaining energy in this period is almost entirely kinetic. This can be understood 
qualitatively by putting cu — —i'j and writing equation (31) in the form 

where primes denote derivatives with respect to r. Assuming that / tends to a 
finite value as r approaches the singular point (7/5)^/^, as in the case of the finite 
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real frequency continuum modes discussed by McClements & Thyagaraja (2004), we 
infer from equation (35) that whereas the magnetic field perturbation /' diverges 
logarithmically in the neighborhood of r = (7/5)^/^, the current (r/')'/r and hence the 
velocity (cf. eq. [16]) diverge as l/[r — (7/S')^/^]. For these continuum eigenf unctions, 
the field energy is thus square integrable whereas the kinetic energy is not. The energy 
components in the late non-oscillatory phase of the resistive MHD simulations can, in 
principle, be synthesized using these singular eigenmodes. In view of the nature of the 
singularities in /' and (r/')'/r, it is not surprising that the kinetic energy, although 
necessarily finite (as indicated by eq. [21], the total energy can never increase), is large 
compared to the field energy. 

3.4 Finite CoUisionless Skin Depth Solutions 

Wc have obtained numerical solutions of equations (15) and (16) with 5^ = 0.01 and S 
ranging from 10^ to 10^, using both types of initial condition discussed in the previous 
two subsections. Figure 6 shows the evolution of the field energy, ion kinetic energy 
and electron kinetic energy in each case, with the left hand and right hand plots 
corresponding respectively to the spatially extended initial condition (eq. [25]) and the 
strongly localized initial condition (eq. [28]). Comparing Figure 6 with Figures 3 and 
4, it can be seen that electron inertia has a negligible impact on the evolution of £f, 
£]^i and £ whenever SS^ < 1: the black, red and blue curves in Figure 3(b) and the 
left hand frame of Figure 6(a), for example, which were obtained using the same initial 
conditions and S = 10^, are indistinguishable. This is particularly noteworthy since 
the results shown in Figure 6 and those discussed in §§3.3 were obtained using two 
independent codes (see Appendix B). When 36^ = 0.01, £ke remains below the other 
energy components at all times for both sets of initial conditions (Fig. 6 [a]). When the 
initial field perturbation has the profile sinTrr, the ratio £ke/£f is well- approximated 
at t = by the equals sign in equation (A7). The ratio is somewhat higher in the 
case of the more localized initial field perturbation: this appears to be due to a greater 
contribution of eigenvalues > Aq in equation ( A8) . Towards the end of simulations 
with SSl — 0.01, £ke becomes comparable to the field energy, but the system evolution 
is still well- approximated by resistive MHD. 

The situation is completely different when S5l > 1: this regime is represented by 
Figures 6(d) and 6(e), with Figure 6(c) representing a transitional case {S5l — 1). 
The oscillatory phase prevalent in the resistive MHD simulations is still present when 
36^ = 1 but wholly absent when this quantity is significantly greater than unity. After 
a few Alfven times, the total energy is dominated by ions and electrons, with S equally 
partitioned between the two species. The total energy still decays, as it must according 
to equation (21), but at a diminishing rate as S increases. In fact, the asymptotic decay 
time of the total energy is equal to 36'^. Given that S is dominated by electrons and 
ions, with £ki — £ke so that £ke — ^/2, the decay rate follows very simply from equation 
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(21): 

McClements & Thyagaraja (2004) demonstrated that the discrete X-point spectrum 
ceases to exist when the intrinsic damping rate of the finite frequency continuum is less 
than that of the least damped discrete mode. As noted previously, the latter varies 
as l/(lnS')^ and is of order unity in Alfven units for low S. Thus, when SS'^ ^ 1 the 
spectrum is purely continuous. The decay rate indicated by equation (36) is exactly 
twice the intrinsic damping rate of finite real frequency continuum eigenmodes (cf. eq. 
[3]). This is to be expected if the energy of the system is assumed to lie predominately 
in the finite frequency continuum, since the ion and electron kinetic energies scale 
as the squares of, respectively, the fluid velocity and the current. As noted in §1, the 
energy decay time SS^ is simply the electron collision time if S is determined by Spitzer 
resistivity. 

For these simulations there is a period in which the field energy decays much more 
rapidly than this, typically on an Alfven timescale: this causes to be energetically 
insignificant at later times. As in the later stages of the resistive MHD simulations, 
the dominance of kinetic energy over field energy when SS"^ > 1 is associated with the 
fact that for continuum eigenfunctions the singularity in the field is square integrable 
while those in the current and velocity are not (McClements & Thyagaraja 2004). The 
only difference in this case is that part of the continuum has finite real frequency. 

The field energy decays on the Alfven timescale even in the collisionless limit {SS"^ — > 
oo): we thus observe irreversible behavior, despite the fact that total energy is 
conserved in this case. As demonstrated by McClements & Thyagaraja (2004), the 
spectrum in the limit SS^ — > cxd consists of a band-limited continuum. The decay 
of field energy in this limit can be interpreted as continuum damping. It is exactly 
analogous to Landau damping: the Vlasov-Poisson system of equations, like the system 
represented by equations (15) and (16) in the limit SS^ — > oo, is conservative but 
exhibits irreversible behavior (see e.g. van Kampen & Felderhof 1967). The analysis 
of inviscid plane Couette fiow by Case (1960) shows that this type of behavior is 
not restricted to plasma kinetic theory or wave-particle interactions but is a general 
consequence of phase mixing. Thyagaraja et al. (2002) have demonstrated that even 
the simplest advection-diffusion equation in one space dimension can exhibit such 
dissipationless damping. In the solar context, phase mixing of MHD waves has been 
investigated by many authors as a possible coronal heating mechanism (e.g. Heyvaerts 
& Priest 1983; Voitenko & Goosens 2000; Tsiklauri, Nakariakov, & Rowlands 2003). 
Porcelli et al. (2002) have shown that the channeling of field energy into kinetic energy 
through phase mixing occurs in the collisionless evolution of magnetic islands. In the 
case of finite S with SS^ > 1, as in Figures 6(d) and 6(e), the field energy decays 
through a combination of phase mixing and the intrinsic finite frequency continuum 
mode damping given by equation (3). 
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Another characteristic of the S5l > 1 regime indicated by Figure 6 is that the evolution 
of S is relatively insensitive to the initial conditions: in Figures 6(d) and 6(e), with SS"^ 
equal to 10 and 100 respectively, the left hand plots are almost identical to the right 
hand plots. In the case of the more localized perturbation, the electrons and ions take 
somewhat less time to reach cquipartition, but the subsequent evolution of the energy 
components is essentially independent of the initial field profile. In the resistive MHD 
regime, in contrast, it is clear from Figures 3 and 4 that the total decrement in the 
energy after a given time interval depends critically on the initial conditions, although 
the rate of energy loss does not. 

The profiles of dt/j/dr and v also depend critically on whether S6'^ is less than or greater 
than unity. Figure 7 shows the evolving velocity profile up to t = 4 for the case of 
initial V oc - exp[-(lnr) VAm^], 5^ = 0.01 and (a) S = 10^ (b) S = 10^ For 35^ < 1, 
the profile steepening observed in the ideal phase (Figs. [1] and [2]) ceases due to 
resistive effects, and the field and velocity are damped without the shapes of their 
spatial profiles being radically altered. A low amplitude ripple appears in the wake of 
the inward propagating pulse, apparently due to finite Sg (the ripple is absent in the 
resistive MHD limit), but, as indicated by Figure 8(a), the overall profile shape is similar 
to the ideal MHD profiles shown in the right hand frame of Figure 2. When 36^ > 1, 
on the other hand, the field and velocity profiles are characterized by short wavelength 
features that persist for long times (Figs. 7[b] and 8[b]). Porcelli et al. (2002) observed 
a similar filamentation process in coUisionless simulations of a magnetic island. The 
appearance of progressively more nodes in the field profile means that the Raylcigh 
quotient defined by equation (A2) increases with time: this is linked to the rapid decay 
of field energy relative to electron kinetic energy noted above. Figures 7(b) and 8(b) 
indicate that the radial wave number spectrum is cascading into progressively shorter 
wavelengths. Examples of such direct cascades investigated previously in the plasma 
literature (e.g. by Thyagaraja, Loureiro, & Knight 2002) show that they play a crucial 
role in the temporal evolution of field energy. 

The occurrence and persistence of small spatial scales for 36"^ > 1 is linked to the 
fact that the only eigenmodes in this limit are singular, and are thus characterized by 
arbitrarily short wavelengths: the energy in the initial perturbation is channeled into 
these short wavelengths, giving rise to the observed profiles. McClements & Thyagaraja 
(2004) showed that the finite frequency continuum is present even when 35^ < 1, but 
in this case the component of the energy in this continuum rapidly damps out (on 
the 36'^ timescale), and the remaining energy decays on a timescale determined by the 
discrete mode damping rate. Since the discrete modes are non-singular, they are not 
characterised by short wavelengths. Thus, in the resistive MHD regime we observe 
relatively smooth spatial profiles (Figs. 7 [a] and 8 [a]) but highly structured evolution 
of the total energy (Figs. 3, 4, 6 [a] and 6[b]), whereas in the 36'^ > 1 regime the 
profiles have fine scale structure (Figs. 7[b] and 8[b]) but the temporal behavior of 
the total energy is relatively smooth (Figs. 6[d] and 6[e]). In the latter case short 
wavelength structures also appear in the current profile and the longitudinal component 
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of the electric field, —dip/dt. Electrons (and ions) would then be accelerated in 
nested cylindrical shells, centered on the X-line, with acceleration occurring in opposite 
directions in adjacent shells. One attractive feature of this scenario is that it suggests 
the possibility of producing large numbers of hard X-ray emitting electrons without 
the creation of large net beam currents, and hence unacceptably large magnetic fields 
(see e.g. Holman 1985). 

McClements & Thyagaraja (2004) noted that the discrete X-point spectrum found 
by Craig & McClymont (1991) is replaced in the ideal limit with a continuum: the 
eigenmodes of this continuum, like those of the non-ideal continua discussed above, are 
singular. The cascading of field and kinetic energy into small spatial scales apparent in 
Figures 1 and 2 is associated with the singular nature of the ideal MHD eigenmodes. 
When S is finite but 5^ = 0, the finite frequency ideal MHD spectrum becomes discrete 
and the eigenmodes non-singular. However, as discussed in §§3.3, the system evolution 
at high S is still accurately described by ideal MHD for several Alfven periods after 
t — 0: energy cascades into small scales in this period, despite the resolution of the 
continuum into discrete modes. Similarly, although one would expect the non-ideal 
continua that exist for finite S and Se to be resolved into discrete spectra by, for 
example, finite gyro radius effects, such effects will not necessarily prevent filamentation 
or the decay of field energy relative to kinetic energy. On the other hand, gyro radius 
effects are likely to impose a lower hmit on the length scale of the filaments: in the 
simulations reported here with SS"^ > 1, there appears to be no such lower limit. 

4 Summary and Application to Solar Flares 

We have studied the relaxation of non-potential perturbations to a current-free 
magnetic X-point, taking into account the effects of resistivity and electron inertia. The 
latter has been shown to have a negligible effect on the evolution of the system whenever 
the collisionless skin depth is less than the resistive scale length. Non-potential 
magnetic field energy in this resistive MHD limit initially reaches equipartition with 
flow energy, in accordance with analytical results obtained using ideal MHD, and is then 
dissipated extremely rapidly, on an Alfvenic timescale that is essentially independent 
of Lundquist number, with the energy being attenuated to a greater extent if the 
initial perturbation is highly localized. Following this period of rapid dissipation in the 
resistive case, the magnetic field energy and kinetic energy decay on a longer timescale 
and exhibit oscillatory behavior: Craig & Watson (1992), who also observed such 
behaviour in resistive MHD X-point simulations, noted that it arises from the existence 
of discrete normal modes with finite real frequency. Eventually, the oscillations 
disappear, and the energy continues to decay on a yet longer timescale. This non- 
oscillatory decay arises from the fact that the discrete modes observed in the earlier 
evolution do not constitute a complete set, even in the resistive limit; the spectrum 
also contains a purely damped continuum. When the collisionless skin depth exceeds 
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the resistive scale length, the system again evolves initially according to ideal MHD. At 
the end of this ideal phase, the field energy decays typically on an Alfvenic timescale, 
while the kinetic energy (which is equally partitioned between ions and electrons in 
this case) is dissipated on the electron coUision timescale. The oscillatory decay in the 
energy observed in the resistive case is absent, but short wavelength structures appear 
in the field and velocity profiles, suggesting the possibility of particle acceleration in 
oppositely-directed current channels. 

We now consider the possible application of these results to the problem of short 
timescale energy release in solar fiares. For a system size Rq equal to the typical 
length of a flaring coronal loop (~ lO^m), Bq ~ O.IT, n = lO^^m-^ T ~ 3 x lO^K, 
and Spitzer resistivity, the Lundquist number S is of the order of 3 x 10^^ while the 
dimcnsionless collisionless skin depth is about 10^^ and hence SS"^, ~ 1. This indicates 
that the collisionless skin depth 6^, although very small compared to the typical spatial 
extent of a solar flare, is actually comparable to the resistive length scale l/S*^/^ if the 
latter is assumed to be determined by Coulomb collisions. Since Spitzer resistivity is 
often invoked as the source of magnetic held energy dissipation in flares, the figures 
quoted above suggest that electron inertial effects are likely to play a role in this 
process. However, in the absence of precise knowledge of magnetic field structures in 
the pre-flare corona, there is no particular reason for identifying Rq with a macroscopic 
flare dimension, given that no other length scales appear in the current-free X-point 
fleld deflned by equation (7). Indeed, it is unlikely that this equation would be a good 
approximation to the coronal field over length scales as large as lO^m. More realistically, 
it could be regarded as a model of a spatially restricted region, with dimensions much 
smaller than a typical fiare size, for which the approximation of a current-free two- 
dimensional X-point is appropriate. Non-potential magnetic fleld energy originating 
from a much more extensive region could still be channeled into the X-point through 
the focusing effect noted by Craig & Watson (1992). A key point to note here is 
that S is defined in terms of a field component which increases linearly with distance 
from a two-dimensional null (it is independent of any longitudinal field component). 
If the boundary is chosen to be arbitrarily close to the origin, the Lundquist number 
is then arbitrarily small. In this specific context, it is thus entirely appropriate to 
consider values of S much smaller than the global Lundquist numbers characteristic of 
the flaring solar corona, and values of Sg much larger than the flgure of 10~^ quoted 
above. 

The choice of Rq is not completely arbitrary, since the boundary conditions imposed at 
this radius must have some effect on the evolution of £. However, the results presented 
in the previous section and by Craig & Watson (1992) suggest strongly that any such 
effect is very weak. In the resistive regime, our results are essentially identical to those 
obtained by Craig & Watson despite the use of different boundary conditions. This 
is consistent with the fact that the time-evolving field and velocity profiles in both 
the resistive MHD and electron inertial regimes invariably show dissipation occurring 
much closer to the X-point than the boundary at R — Rq. 
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Regardless of the initial conditions, our results show the greater part of the non- 
potential magnetic field energy either being dissipated or, in the case of S5l > 1, 
being converted into kinetic energy, on the Alfven time ta — Rq/caq- In the event 
of S5l being greater than unity, with S determined by Spitzer resistivity, the kinetic 
energy is dissipated on the electron collision time Xg. Since S5l is simply the ratio 
Te/TA (cf. eqs. [3] and [4]), the overall energy dissipation timescale can thus be written 
as max(rA,Te). As noted in §1, is typically less than a second in the flaring corona. 
The Alfven time could also be less than a second, if the X-point field were sufficiently 
highly sheared. Again putting n — lO^^m"^, and assuming that the field rises to, say, 
O.OIT at a radial distance from the null of lO^m, we obtain ta — 0.1s. Since it is the 
destruction of magnetic flux that gives rise to an accelerating electric field in the z- 
direction = —dip/dt, the key timescale from the point of view of energetic particle 
production is always ta, provided that the acceleration time is no longer than this. 
The latter timescale depends, of course, on the amplitude of the initial magnetic field 
perturbation. Since equations (15) and (16) are linear, the perturbation amplitude is 
arbitrary as far as the results shown in Figures 1-8 are concerned. However, we can at 
least check that the quantity of magnetic flux A^^ that must be destroyed to produce, 
say, a 30keV electron in 0.1s is broadly consistent with the linear approximation. The 
energy £e acquired by a nonrelativistic electron in time t is 

4 - V-^. (37) 

so that 

AV^ = (2m£e)'/'^. (38) 

The perturbed field is ~ /^-^/Rq and the unperturbed field is Be ^ Bq: their ratio 
is 

Using the parameter values invoked above, we obtain B^/Be ~ 10~^. Even allowing 
for steepening of the perturbed field during the course of each simulation, this crude 
estimate suggests that the production of SOkeV electrons on sub-second timescales is 
easily compatible with the use of a linear model. 

We have omitted from our model many effects that could in principle infiuence the 
field decay time. These include ion gyro radius effects, equilibrium currents and fiows, 
longitudinal magnetic fields, Hall currents, nonlinear effects, pressure gradients, the 
effects of three dimensional geometry, compressibility and kinetic effects. For example, 
both McClymont & Craig (1996) and Biskamp et al. (1997), using very different 
physical models, have found that the reconnection rate can be strongly affected by 
the presence of a longitudinal magnetic field. Our purpose in this paper has been 
to achieve a reasonably complete understanding of a relatively simple physical model, 
with the minimum number of free parameters. To achieve this understanding, we have 
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made full use of the spectral analysis presented in our earlier paper (McClements & 
Thyagaraja 2004), although the complementary initial value approach was essential in 
order to demonstrate some of our key results, notably coUisionless continuum damping 
of the field energy on the Alfven timescale. The scheme we have developed could be 
readily extended to incorporate most of the additional effects listed above. The results 
in this paper demonstrate that resistivity and electron inertia alone are sufficient to 
endow perturbed magnetic X-points with a rich physical structure. 

This work was funded by the United Kingdom Engineering and Physical Sciences 
Research Council. Helpful conversations with Per Helander are gratefully 
acknowledged. 



Appendix A: Relation between Electron Kinetic 
Energy and Field Energy 



We can establish a useful relation between the field energy 8i and the electron kinetic 
energy as follows. The latter is given by 
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where </? = —dip/dr is the (purely azimuthal) perturbation to the magnetic field. This 
expression prompts us to consider the Rayleigh quotient R[(p\ defined by the expression 
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Following the ideas underlying Wirtinger's inequality (Hardy, Littlewood, & Polya 
1967), we attempt to minimize this functional over all test functions (p that vanish at 
r — and r — 1. According to the Rayleigh- Ritz principle frequently employed in wave 
mechanics (e.g. Mathews & Walker 1964), the minimum is attained for (p satisfying 
the Euler-Lagrange equation 
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where Aq > is the smallest eigenvalue of this equation. Denoting the corresponding 
eigenfunction by (fimm, it is straightforward to establish that 



AJi(Aor), 



(A4) 



where Ji is the Bessel function of order one and A is a normalization constant that 
can be chosen to be unity without loss of generality. Evidently Aq is equal to the 
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first positive zero of Ji, i.e. ji i ~ 3.8 (e.g. Abramowitz & Stegun 1965). From 
equations (A2) and (A3), integrating the numerator of the former by parts and using 
the boundary conditions (p — ai r = and r = 1, we infer the inequahty 

R[^] >Xl^ Jti - 14-4. (A5) 

The equals sign in this result applies if and only ii (p — (puiin — Jii^or). Since the field 
energy is given by 

^ f dip\^ 1 



^'=2 Jo [d^^ 
it follows that R[(p\ — E^ej (^f(^e) hence 

4e > \l^lEi. (A7) 

Thus, for any perturbation to the equilibrium X-point magnetic field, the electrons 
must have finite kinetic energy at all times. 

In general, </? can be expanded in a Fourier- Bessel series of the form 

oo 

ip=Y. x/2a„Ji(A„r)/J2(A„), (A8) 

n=0 

where A„ is the (n + l)-th positive zero of Ji, the a„ are constants, and we have used 
the relation (Whittaker & Watson 1927) 
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Ji{Xnr)Ji{Xmr)rdr = J^{Xn)6mn/'2, (A9) 



J2 being the Bessel function of order two and the Kronecker delta symbol, to ensure 
that the eigenf unctions form a complete orthonormal set. The Rayleigh quotient is then 
given by 



Y^OD 2 \2 



m = " (Aio) 



i.e. a weighted mean eigenvalue of the differential operator in equation (A3), the 
weights in question being the squares of the Fourier- Bessel amplitudes a„ of the 
perturbed field (p. Since the A„ form a monotonic increasing sequence and the index n 
counts the number of nodes of the eigenfunction, it follows that R can be regarded as 
a measure of the effective radial wave number of p (cf. Thyagaraja 1979). The initial 
value of this measure relates to the localization of the initial state. If (p became more 
localized or developed more nodes as the system evolved, it would be characterized 
by relatively high eigenvalues A^, and we would then expect a relatively high ratio of 
electron kinetic energy to field energy. 
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Appendix B: Numerical Algorithms 



To solve equations (15) and (16) numerically it is convenient to introduce new variables 
u — rv, b — rdip/dr and 
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Differentiating equation (15) with respect to r and multiplying by r we obtain 

dw du 1 d ( 1 dh\ 
dt dr S dr \r dr j ^ 
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while the momentum equation becomes 
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For the case of finite 6e we solve these equations using a staggered leap-frog method, 
in which u{r,t) is displaced with respect to b{r,t) and w{r,t) by half the space step, 
Ar. Labelling the time step by At, and the space and time grid points by i and n 
respectively, we approximate equations (Bl) - (B3) using the finite difference scheme 
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Equations (B4) - (B6) are advanced in time using an iterative predictor /corrector 
scheme, with a tri-diagonal matrix algorithm being used to solve equation (B6) for b 
in terms of w. 

The resistive solutions {5e = 0) presented in §§3.3 were obtained using a code based 
on a somewhat simpler technique. The discrete values of the field variable b are 
again staggered in space with respect to those of the velocity variable u, but the 
scheme is explicit, and thus does not require iteration or matrix inversion. As noted in 
§§3.4, results obtained using the two codes are essentially identical whenever SS"^ < 1. 
Moreover, for large S and early times, both schemes recover the ideal solutions given 
by equations (26), (27) and (29). 
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Figure 1: Left plot: time evolution of field energy (red curve) and ion kinetic energy 
(blue curve) in ideal MHD limit for dip /Or oc sin(7rr), f = at t = 0. Right plot: time 
evolution of dip /dr. The curves correspond to t = (black), t = 1 (red), t = 2 (green), 
and t = 3 (blue). 
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Figure 2: Left plot: time evolution of field energy (red curve) and ion kinetic energy 
(blue curve) in ideal MHD limit for ip oc — exp[— lOO(lnr)^], f = at t = 0. Right 
plot: time evolution of dip /dr. The curves correspond to t = (black), t = 1 (red), 
t = 2 (green), and t = 3 (blue). 
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Figure 3: Time evolution of normalized total energy S/Sq (black curves), field energy 
(red) and kinetic energy (blue) for 4 = and ^ = 10 (a), 10^ (b), 10^ (c), 10^ (d), 10^ 
(e), 10^ (f), 10'' (g), 10*^ (h). The initial conditions are v = and dip/dr oc sin(7rr) 
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Figure 5: Time evolution over 200 Alfven periods of normalized total energy EjE^ 
(black curves), field energy (red) and kinetic energy (blue) for 5e = and S = 10*^ (a), 
10^ (b). The initial conditions are v = and dip /Or oc sin(7rr) 



27 



o 



(a) 



dip/dr oc sin(7rr) 
at t = 




\ -0 oc 


-exp[-100(lnr)2] ; 




at t = ; 






^ vv — ^ 

; (a) 






(d) 



(e) 



Figure 6: Time evolution of normalized total energy S/So (black curves), field energy 
(red), ion kinetic energy (blue) and electron kinetic energy (green) for Se = 0.01 and 
(a) S = 10^, (b) 10^, (c) 10^ (d) 10^ (e) 10^. The initial conditions are w = and 
dip/dr oc sin(7rr) (left plots), ip — exp[— lOO(lnr)^] (right plots). 
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Figure 8: Velocity profile at t = 4 in simulations with 5e = 0.01 and (a) S = 10^, (b) 
S — 10^. The initial conditions are v — and t/j oc — exp[— lOO(lnr)^]. 
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